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CHAPTER SEVEN 
INTERNAL WAVES IN A STRATIFIED FLUID 

1 Introduction. 

The atmosphere and ocean are continuously stratified due to change in temperature, 
composition and pressure. These changes in the ocean and atmosphere can lead to 
significant variations of density of the fluid in the vertical direction. As an example, 
fresh water from rivers can rest on top of sea water, and due to the small diffusivity, the 
density contrast remains for a long time. The density stratification allows oscillation of 
the fluid to happen. The restoring force that produces the oscillation is the buoyancy 
force. The wave phenomena associated with these oscillations are called internal waves 
and are discussed in this chapter. 

2 Governing Equations for Incompressible Density- 
stratified Fluid. 

We are going to derive the system of equations governing wave motion of an incompress- 
ible fluid with continuous density stratification. Cartesian coordinates x, y and z will 
be used, with z measured vertically upward. The velocity components in the directions 
of increasing x, y and z will be denoted as u, v and w. The fluid particle has to satisfy 
the continuity equation 

1 Dp du dv dw 
p Dt dx dy dz 

and the momentum equations 
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du du dp 

^ dt dx dx ' 
dv dv dp 

^ dt dy dy ' 

dw dw dp 

P ^ + W ^ = -d- Z -^ (2 - 4) 

where p and p are, respectively, the fluid density and pressure. The fluid is taken to 
be such that the density depends only on entropy and on composition, i.e., p depends 
only on the potential temperature 9 and on the concentrations of constituents, e.g., the 
salinity s or humidity q. Then for fixed 9 and q (or s), p is independent of pressure: 



p = p(9,q). (2.5) 

The motion that takes place is assumed to be isentropic and without change of phase, 
so that 9 and q are constant for a material element. Therefore 

Dp = dpD9 dpDq = 

Dt d9 Dt dqDt [ ' ' 

In other words, p is constant for a material element because 9 and q are, and p depends 
only on 9 and q. Such a fluid is said to be incompressible, and because of (2.6) the 
continuity equation (2.1) becomes 

du dv dw 
dx dy dz 

For an incompressible fluid, the density p satisfies the density equation 

1^ = 0. (2.8) 
pDt v ' 

Assuming that the velocities are small, we can linearize the momentum equations to 
obtain 
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9u dp 

dv dp , 

»m = ~iy (2 ' 10 » 

dw dp . . x 

Next, we consider that the wave motion results from the perturbation of a state of 
equilibrium, which is the state of rest. So the distribution of density and pressure is the 
hydrostatic equilibrium distribution given by 



I = -»* < 2 - 12 > 

When the motion develops, the pressure and density changes to 



p=p(z)+p', (2.13) 
p=p(z)+p', (2.14) 

where p' and p' are, respectively, the pressure and density perturbations of the "back- 
ground" state in which the density p and the pressure p are in hydrostatic balance. The 
density equation now assume the form 



f + u 9 / + v a / + w f + w a / = . ( 2 . 1B ) 

dt dx dy dz dz 

The nonlinear terms u(dp'/dx), v(dp'/dy) and w(dp'/dz) are negligible for small am- 
plitude motion, so the equation (2.15) simplifies to 



which states that the density perturbation at a point is generated by a vertical advection 
of the background density distribution. The continuity equation (2.7) for incompressible 
fluid stays the same, but the momentum equations (2.9) to (2.11) assume the form 
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p m = "to* (2 - 17) 

<"»> 

We would like to reduce the systems of equations (2.7), (2.16) and (2.17) to (2.19) to a 
single partial differential equation. This can be achieved as follows. First, we take the 
time derivative of the continuity equation to obtain 



d 2 u d 2 v d 2 w „ _ 

+ TTTTj + 77777" = 0. (2.20) 



. d 2 u 


d 2 p' 


dxdt 


dx 2 ' 


_d 2 v 


d 2 P ' 


dydt 


dy 2 ' 


d 2 w 


d 2 p' 


P dt 2 


dtdz 



dtdx dtdy dtdz 

Second, we take the x,y and t derivatives, respectively, of the equations (2.17) to (2.19) 
and we obtain 



p dxm = " a^' (2 - 21) 

d 2 v d 2 p' 

P—, = ( 2 - 22 ) 

<"»> 

If we substitute equations (2.21) and (2.22) into equation (2.20), we obtain 
We can eliminate p' from (2.23) by using equation (2.16) to obtain 



d 2 w d 2 p' dp 

*w = -aifc + >&°- <2 ' 25) 

Third, we apply the operator + to equation (2.25) to obtain 

_d 2 (d 2 w d 2 w\ d 2 (d 2 p' d 2 p'\ dp(d 2 w d 2 w\ 

P W\^ + ^) = ~^z\^^^) +9 Yz\^ + ~d^)' ' 
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Next, we use equation (2.24) to eliminate p' from equation (2.26), which gives the 
following partial differential equation for w: 



dt 2 



d 2 w d 2 w 



dx 2 



dy 2 



18_ 

p dz 



_dw 
dz 



N 2 



d 2 w d 2 w 



dx 2 



dy 2 



(2.27) 



where we define 



N 2 (z) = - y -^, (2.28) 
pdz 

which has the units of frequency (rad/sec) and is called the Brunt- Vaisala frequency or 
buoyancy frequency. If we assume that w varies with z much more rapidly than p(z), 
then 



pdz \ dz J dz 2 ' 



and (2.27) can be approximated by the equation 



d 2 ( d 2 w d 2 w d 2 w\ , t9 / d 2 w d 2 w\ 

5?(a? + W + aJ) +N vj = °' <2 ' 30) 

The assumption above is equivalent to the Boussinesq approximation, which applies 
when the motion has vertical scale small compared with the scale of the background 
density. It consists in taking the density to be constant in computing rates of change 
of momentum from accelerations, but taking full account of the density variations when 
they give rise to buoyancy forces, i.e., when there is a multiplying factor g in the ver- 
tical component of the momentum equations. The Boussinesq approximation leads to 
equation (2.30) for the vertical velocity w. 



3 The Buoyancy Frequency (Brunt- Vaisala frequency). 

Consider a calm stratified fluid with a static density distribution p(z) which decreases 
with height z. If a fluid parcel is moved from the level z upward to z + (, it is surrounded 
by lighter fluid of density p(z + £). The upward buoyancy force per unit volume is 
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g\p{z + Q-p{z)]ngfy, (3.31) 
and it is negative. Applying Newton's law to the fluid parcel of unit volume, we have 



9 2 C dp > 

or 



dt 



^ + N 2 ( = 0, (3.33) 



where 



N 2 {z) = _gtf (334) 

p dz 

which is called the buoyancy frequency or the Brunt Vasala frequency. This elementary 
consideration shows that once a fluid is displaced from its equilibrium position, gravity 
and density gradient provide restoring force to enable oscillations. 

4 Internal Gravity Waves in Unbounded Stratified 
Fluid. 

Consider the case in which the buoyancy (Brunt-Vaisala) frequency N is constant 
throughout the fluid. Traveling wave solutions of (2.30) can be found of the form 



w = wo cos(kx + ly + mz — cot), (4.35) 

where w is the vertical velocity amplitude and it = (k, I, m) is the wavenumber of the 
disturbance, and u is the frequency. In order for (4.35) to satisfy the governing equation 
(2.30) for the vertical perturbation velocity, u and Jz must be related by the dispersion 
relation 
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2 (P + / 2 )iV 2 

w = 75 79 9 - 4.36 

k 2 + I 2 + m 2 v ' 

Thus internal waves can have any frequency between zero and a maximum value of 
N. The dispersion relation for internal waves is of quite a different character compared 
to that for surface waves. In particular, the frequency of surface waves depends only 
on the magnitude of the wavenumber, whereas the frequency of internal waves is 
independent of the magnitude of the wavenumber and depends only on the angle <f> that 
the wavenumber vector makes with the horizontal. To illustrate this, we consider the 
spherical system of coordinates in the wavenumber space, namely, 



k =|T| cos(^) cos(6>) (4.37) 

/ =|~^l C0S 0) sm (#) (4.38) 
m=\~t\sm(cf)) (4.39) 

The coordinate system in the wavenumber space is given in the figure 1. 
The dispersion relation given by equation (4.36) reduces to 



u 2 = Ncos(4>). (4.40) 

Now we can write expressions for the quantities p', p', u and v. From equation (2.20) we 
can write 



1 fd 2 p' d 2 p'\ d 2 w 

+ — — - I = = urniwo cos(kx + ly + mz — cot), 



po \ dx 2 dy 2 J dtdz 
which implies that the perturbation pressure p' is given by 



. umw po ,aa-,\ 
P = - ^ 2 - py /2 cos(kx + ly + mz - cot). (4-41) 

From equation (2.16) we have the perturbation density p' given by 
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Figure 1: Coordinate system in the wavenumber space. 
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/7V 2 \ 

p' = — — PoWq sin(kx + ly + — wt). (4.42) 

The horizontal velocity components can be found from equations (2.17) and (2.18), 
which give 



{u, v) = — {k, l)(k 2 + I 2 ) 1 mw cos(kx + ly + mz — cot) (4.43) 
= (M)(wp )-y. (4.44) 

The above relations between pressure and velocity fluctuations can be useful for de- 
ducing wave properties from observations at a fixed point. For instance, if the horizontal 
velocity components and perturbation pressure of a progressive wave are measured, the 
horizontal component of the wavenumber vector can be deduced from (4.44). 

A sketch showing the properties of a plane progressive internal wave in the vertical 
plane that contains the wavenumber vector is presented in figure 2. The particle motion 
is along wave crests, and there is no pressure gradient in this direction. The restoring 
force on a particle is therefore due solely to the component g cos <f> of gravity in the 
direction of motion. The restoring force is also proportional to the component of density 
change in this direction, which is cos <f)j^ per unit displacement. 

Consider now the succession of solutions as (f> progressively increases from zero to 
7r/2. When 4> = 0, a vertical line of particles moves together like a rigid rod undergoing 
longitudinal vibrations. When the line of particles is displaced from its equilibrium, 
buoyancy restoring forces come into play just as if the line of particles were on a spring, 
resulting in oscillations of frequency N. The solution for increasing values of 4> corre- 
spond to lines of particles moving together at angle <f> to the vertical. The restoring force 
per unit displacement (cos 4>dp' /dz) is less than the case where 4> = 0, so the frequency 
of vibration is less. As <f> tends to tt/2, the frequency of vibration tends to zero. The 
case 4> = 7r/2 is not an internal wave, but it represents an important form of motion that 
is often observed. For instance, it is quite common on airplane journeys to see thick 
layers of cloud that are remarkably flat and extensive. Each cloud layer is moving in its 
own horizontal plane, but different layers are moving relative to each other. 
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Figure 2: The instantaneous distribution of velocity, pressure, and buoyancy perturba- 
tions in an internal gravity wave. This is a view in the x, z plane. The phase of the wave 
is constant along the slanting, dashed, and solid lines. Velocity and pressure perturba- 
tions have extrema along the solid lines; buoyancy perturbations are zero along the solid 
lines. Buoyancy perturbations have extrema, and velocity and pressure perturbations 
are zero along dashed lines. Small arrows indicate the perturbation velocities, which are 
always parallel to the lines of constant phase. Large heavy arrows indicate the direction 
of phase propagation and group velocity. 
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4.1 Dispersion Effects. 

In practice, internal gravity waves never have the form of the exact plane wave given 
by equation (4.35), so it is necessary to consider superposition of such waves. As a 
consequence, dispersion effects become evident, since waves with different frequencies 
have different phase and group velocities as we are going to show in this section. For 
internal waves, surfaces of constant frequency in the wavenumber space are the cones 
4> = constant. The phase velocity is parallel to the wavenumber vector and it lies on a 
cone of constant phase. Its magnitude is 




(4.45) 



The group velocity C g is the gradient of the frequency u in the wavenumber space and 
therefore is normal to the surface of constant frequency u. It follows that the group 
velocity is at right angles to the wavenumber vector. When the group velocity has an 
upward component, therefore, the phase velocity has a downward component, and vice 
versa. The group velocity vector is 



C g = -=p sin </»(sin 0cos 9, sin 4> sin 9, — cos </>). (4.46) 
| k | 

Therefore, the magnitude of the group velocity is (y^j) sin</>, and its direction is at 
an angle <f> to the vertical, as illustrated in the figure 3. 

To illustrate the effects of dispersion, we consider the case of two dimensional mo- 
tions. We consider only the coordinates x and z. In this case, the wavenumber is the 
vector (k,m). We consider an initially localized wave packet. Due to dispersion effects, 
the wave packet spreads and moves according to the group velocity vector C g , which 
now simplifies to 



C g = -r=y- sin </>(sin </>, — cos </>). (4.47) 
| A: | 

The phase velocity is perpendicular to the group velocity vector, so the wave crests 
(lines of constant phase) move perpendicularly to the direction of propagation of the 
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Figure 3: Wavenumber vector and group velocity vector. 
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wave packet. The phase velocity is given by the equation (4.45), where the wavenumber 
vector Jz makes an angle <f> with the horizontal direction (see figure 1, but now set 
6 = 0). 

To illustrate the effects of dispersion, we consider three different animations of a 
localized wave packet for the density perturbation p' . The perturbation density p' is 
related to the vertical velocity w by the equation 

- = —w, (4.48) 

and the governing equation for the vertical velocity w is given by the equation (2.30). 
To obtain the evolution in time of an initially localized wave packet for the perturbation 
density, we apply a two-dimensional Fourier transform to equations (2.30) and (4.48). 
The two-dimensional Fourier transform pair considered is 

/oo poo 
dx / dz {exp(— ikx — imz)u(x, z)} (4.49) 
oo J — oo 

and 



poo poo 

u(x, z) = / dk I dm {exp(—ikx — imz)u(k,m)} . (4.50) 

^ J —oo J —oo 

The Fourier transform of the equation (2.30) is given by the equation 



dt 2 k 2 + m 2 

which has solution of the form 



8 2 w N 2 k 2 . n fAK ,. 



w(k, m, t) = A(k, m) exp(zwt) + B(k, m) exp(— iut), 
where uj is given by the dispersion relation 



(4.52) 
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The Fourier transform of the equation (4.48) is given by the equation 

dp' p N 2 



dt " g 

From equations (4.51) and (4.54) we have that 



-w. (4.54) 



p'(k, m, t) = — ^- r {— iA(k, m) expUut) + iB(k, m) expUut)} , (4.55) 

gu(k,m) 

where the constants A and B are determined from the Fourier transform of the initial 
conditions for p', given by the equations 



p'(x,z,0)=f(x,z), (4.56) 
dp' 

-£{x,z,0)=0, (4.57) 
which implies that the constants A(k,m) and B(k,m) are given by the equations 



The perturbation density p'(x,z,t) is finally given by the equation 



1 roc roo 

p'(x,z,t) = / dk dm < f(k, m) exp(ico(k, m)t) + f(k, m) exp(— ico(k, m)t) > exp( 

87T^ J_ 0O J x L ' J 

(4.60) 

The function f(x,z) and its Fourier transform f(k,m) are given by the equations 



f(x, z) =-!- exp \—\x 2 a 2 — \z 2 t 2 ) cos(kx + m), (4.61) 
Z7T \ 2 2 J 

s„ , if [ l(k-k) 2 l(m-m) 2 \ ( l(k + k) 2 l(m + m) 2 

f(k, m) = < exp — ^ —!■ ^ —- + exp — ^ ^ —- 

2(jt y y 2 a 2 2 t J \ 2 2 T 

(4.62) 
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In the animation which follows, we show the results from the numerical evaluation of 
the inverse Fourier transform in equation (4.60) for a sequence of values of the variable 
t with / given by equation (4.62). 

The first animation has as initial condition a Gaussian wave packet with a = 1/4, 
t = 1/4 and k = m = |. This initial wave packet has a circular shape and splits in 
two parts as time increases. This two parts propagate in opposite directions from each 
other. Since the x and z components of the main wavenumber are equal and positive 
and the wave packet has the same modulation along the x and z directions (a = r), 
the two parts of the initial wave packet travels towards the middle of the second and 
fourth quadrants, as we see in the animation. For the wave packet in the second (fourth) 
quadrant the group velocity vector points away from the origin towards the middle of 
the second (fourth) quadrant, so the phase velocity, which is orthogonal to the group 
velocity, is oriented in the anti-clockwise (clockwise) sense, as we can see from the crests 
movement in the animation. When the two splited parts of the initial wave packet are 
still close, we see some constructive and destructive interference. To see this animation, 
click here. 

The second animation has as initial condition a Gaussian wave packet with a = 1/2, 
r = 1/100 and k = m = f. This initial wave packet has a shape of an elongated ellipse 
in the x direction. In the movie frame, this initial wave packet looks almost without 
variation in the x direction. The wave packet splits in two parts as time increases. 
This two parts propagate in opposite directions from each other, in a way similar to 
the previous example. The interference effect between the two splited wave packet for 
earlier times is more intense than what was observed in the previous example, as we can 
see in the animation. To see it, click here. 

The third animation has as initial condition a Gaussian wave packet with a = 1/2, 
t = 1/20 and k = m = |. This initial wave packet has a shape of an elongated ellipse 
in the x direction. The movie frame shows the whole wave packet, which splits in two 
parts as time increases. This two parts propagate in opposite directions from each other, 
but the group velocity vector has a slightly smaller component in the x direction. This 
is due to the difference of the modulation of the wave packet in the x and z directions, 
as we can see in the animation. To see this animation, click here. 
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4.2 Saint Andrew's Cross. 

Here we discuss the wave pattern for internal waves produced by a localized source on 
a sinusoidal oscillation, like an oscillating cylinder for example, in a fluid with constant 
density gradient (the buoyancy frequency is constant). For sinusoidal internal waves, 
the wave energy flux 7^ = p'nt (the perturbation pressure p' is given by equation (4.41) 
and the components of the velocity vector are given by equations (4.44) and (4.35)) 
averaged over a period is given by the equation 

-> 1 wlNmpQ a A\ (A R1\ 

* = T^ - To rT~ \sm (p cos 9, sin <p sin 9, — cos (p \ , (4.6dJ 

2 k + 1 2 

which is parallel to the group velocity, according to equation (4.46). Therefore, for 
internal waves the energy propagates in the direction of the group velocity, which is 
parallel to the surfaces of constant phase. This fact means that internal waves generated 
by a localized source could never have the familiar appearance of concentric circular 
crests centered on the source, as we see, for example, for gravity surface waves. Instead, 
the crests and other surfaces of constant phase stretch radially outward from the source 
because wave energy travels with the group velocity, which is parallel to surfaces of 
constant phase. 

For a source of definite frequency u < N (less than the buoyancy frequency), those 
surfaces are all at a definite angle 

<j> = cos"V/iV), (4.64) 

to the vertical] therefore, all the wave energy generated in the source region travels at 
that angle to the vertical. Accordingly, it is confined to a double cone with semi-angle 
4>. The direction of the group velocity vector along the double cone is specified by the 
fact that energy has to radiate out from the source. The direction of propagation of the 
lines of constant phase is also specified in terms of the direction of the group velocity 
and by the fact that the phase velocity 



~(f = — __), - {cos 4> cos 9, cos 4> sin 9, sin </>} (4.65) 
I k I 
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is orthogonal to the group velocity, and that 

-=> N 

C + C g = ^-{cos6, sin6, 0} . (4.66) 
| A; | 

Then, given the direction of the group velocity, the orthogonality of the phase and 
group velocity plus the condition (4.66), the direction of the phase velocity is specified. 
If the group velocity has a positive vertical component, the phase velocity has a negative 
vertical component and vice- versa. The two-dimensional case of an oscillating cylinder 
is illustrated in figure 4. 

This unique property of anisotropy has been verified in dramatic experiments by 
Mowbray and Stevenson. By oscillating a long cylinder at various frequencies vertically 
in a stratified fluid, equal phase lines are only found along four beams forming "St. 
Andrew's Cross", see figure 5 for uj/N = 0.7 and uj/N = 0.9. It can be verified that the 
angles are 4> = 45 degrees for uj/N = 0.7, and 4> = 26 degrees for uj/N = 0.9, in close 
accordance with the condition (4.64). 

5 Waveguide behavior. 

In this section we study free wave propagation in a continuously stratified fluid in the 
presence of boundaries, like an ocean or an atmosphere. Attention is restricted to the 
case in which the bottom is flat, but neither the hydrostatic approximation nor long- 
wave approximation will be made. The equilibrium state that is being perturbed is the 
one at rest, so density, and hence buoyancy frequency, is a function only of the vertical 
coordinate z. We start with an ocean, which has an upper boundary. The atmosphere 
is somewhat different from the ocean since it has no definite upper boundary, so a study 
of waves in this situation is made later in this section. 

5.1 The oceanic waveguide 

Since we assume the undisturbed state as the state of rest, fluid properties are constant 
on horizontal surfaces and, furthermore, the boundaries are horizontal. Solutions of the 
perturbation equation (2.27) can be found in the form 




Figure 4: Phase and group velocities. 




Figure 5: St Andrew's Cross in a stratified fluid. In the top figure uj/N = 0.9 and in 
the left bottom figure uj/N = 0.7. 
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w(x, y, z, t) = w(z) exp[i(kx + ly — cot)] 



(5.67) 



The equation for w(z) can be found by substitution of equation (5.67) into the governing 
equation (2.27). We obtain 



ld_ 

p dz 



_dw 
^ dz 



+ 



(N 2 - u 2 ) 



(k 2 + l 2 )w(z) = 



(5.68) 



The boundary conditions for this equation are the bottom condition of no flux across 
it, given by the equation 



w(z) = at z = —H, 



(5.69) 



and at the free-surface we have the linearized condition 



~dt 



pgw(z) at z = 0, 



(5.70) 



where p' is the perturbation pressure. From this equation we can obtain a free-surface 
boundary condition for w(z). We apply the operator to the equation (2.24), and then 
we substitute equation (5.70) into the resulting equation. As a result, we obtain the 
equation 



d 3 



IV 



dt 2 dz 



' d 2 w d 2 w\ 
9 I tt^t + -rr^- ) at z = 



dx 2 



dy l 



(5.71) 



Now, if we substitute equation (5.67) into the equation (5.71), we obtain the free-surface 
boundary condition for w(z), which follows 



d 2 w (N 2 -oj 2 ) n2 l2 , _ , n 

-(k 2 + l 2 )w(z) = at z = 0. 



+ 



(5.72) 



dz 2 w 

To simplify the governing equation for w(z), we make the Boussinesq approximation, 
such that equation (5.68) simplifies to 
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ff + (^V + ^ W=0 , (5.73) 

OZ z UJ Z 

with boundary conditions given by equations (5.72) and (5.69). The two boundary 
(bottom and free-surface) have the effect of confining the wave energy to a region of 
finite extent, so the ocean can be considered as a waveguide that causes the energy to 
propagate horizontally. 

A useful piece of imaginary is to picture internal waves propagating obliquely through 
the ocean, reflections at the upper and lower boundaries ensuing no loss of energy from 
the wave guide, whereas horizontal propagation is uninhibited. 

Next, we obtain the general solution of equation (5.73) under the boundary condi- 
tions (5.72) and (5.69). We first consider the case where uj 2 > N 2 . For this case the 
general solution has the form 



cosh[m(z + H)] , 2 (uj 2 -N 2 ) 

\) m with m 2 = ± 

cosh(milJ u z 



w(z) = L '";lT r 7 with m 2 = KW ~ iV > {k 2 + 1 2 , (5.74) 



which already satisfies the bottom boundary condition. The free-surface boundary con- 
dition (5.72) gives the dispersion relation 



mtanh(m#) = -^{k 2 + I 2 ), (5.75) 

which is similar to the dispersion relation for surface waves. Actually, the solution (5.74) 
is not an internal wave, but a surface gravity wave. To have internal waves, we need 
that uj 2 < N 2 . This is the next case to consider. We consider the general solution of 
equation (5.73), which is given by the equation 



w(z) = sin[m(z + H)} with m 2 = ± — '-(k 2 + / 2 ), (5.76) 

CO 

which already satisfies the bottom boundary condition. If we substitute equation (5.76) 
into the free-surface boundary condition (5.72), we obtain the dispersion relation 
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N 2 - co 2 = gmtanh(mH). (5.77) 

For a given value of the frequency u, this dispersion relation gives a countable set of 
values for the modulus of the horizontal component (k 2 + I 2 ) of the wavenumber, or for 
a given value of the modulus of the horizontal component of the wavenumber, we have 
a countable set of possible value for the frequency u. For to smaller or of the same order 
of the buoyancy frequency N, the rigid lid approximation can be made, i. e., the left 
hand side of equation (5.72) is small compared with the right hand side, so equation 
(5.72) reduces to 

w(z) = at z = 0. (5.78) 
This boundary condition gives a dispersion relation of the form 



sm(mH) = (5.79) 



or 



2 (k 2 + l 2 )N 2 H 2 , x 

" = n^ + (^ + P)ffl ' n = 1 ' 2 ' 3 '-' (5 - 80) 

which is close to the result given by the dispersion relation given by the free-surface 
boundary condition (5.77). The value of m for the case with a free-surface is slightly 
larger than the case with the rigid lid approximation. 

If the ocean is perturbed with a spatial structure of one of the modes (a specific value 
of m for a given u), then the subsequent behavior in time is described by equation (5.67), 
i. e., there is an oscillation with a particular frequency. Such a situation, however, is 
unlikely, so it is necessary to represent the initial structure in space as a superposition 
of modes (for a given u, we have a countable set of values for k 2 + I 2 ). Then each of 
these will behave in time as found above, and so the solution can be constructed at all 
times by taking the appropriate superposition of modes. 
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5.2 Free Waves in a semi-infinite region. 

The atmosphere does not have a definite upper boundary as does the ocean, so solutions 
of equation (5.73) will now be considered for the case of a semi-infinite domain z > 0. 
In this case there are two types of solutions, the first being typified by the case N = 
constant. The only solutions of equation (5.73) that satisfy the condition at the ground 
z = and remain bounded at infinity are sinusoidal, i. e., 



w(z) = sm(mz), (5.81) 

where m has the same expression as the one given in equation (5.76). There is now 
no restriction on m, so according to the functional relation between m and u given in 
equation (5.76), the frequency u can have any value in the range < u < N, i. e., there 
is a continuous spectrum of solutions. Superposition of such solutions can be used to 
solve initial-value problems, and have the from of Fourier integrals. 

When N varies with z, there is another type of solution possible, namely, one that 
satisfies the condition at the ground yet decays as z — > oo. These are waveguide modes, 
and there are, in general, only a finite number possible. A simple example is provided 
by the case in which a region of depth H of uniform large buoyancy frequency N\ 
underlies a semi-infinite region of uniform small buoyancy frequency N 2 . The layer with 
buoyancy frequency Ni has depth H and lies at < z < H and the semi-infinite layer 
with buoyancy frequency N 2 lies at z > H. For < u < N 2 , the solution in both layers 
has the form given by equation (5.81) with m = m\ in the first layer and m = mi in 
the second layer. The wave frequency is constant across the interface of the two layers, 
which gives the relation 



N\ N, 



ml + k 2 + I 2 ml + k 2 + I 2 (5-82) 
between the vertical wavenumbers nil and m 2 . For this case, the spectrum is continuous 
and co can assume any value between and N 2 . This is not true for the case when 
N 2 < oo < Ni, when the frequency u can assume only a finite set of values in the range 
N 2 < uj < N\. In this case, the solution of equation (5.73) for the first layer is given by 
the equation 
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w(z) = sin(mi^) where ml 



(5.83) 



and in the second layer we have solution given by the equation 



w(z) = exp(— 72) where 7 2 



(^ 2 " N 2 ) 



(k 2 + l 2 ). 



(5.84) 



At the intersection z = H between the two layers, the perturbation pressure p' and the 
vertical velocity w should be continuous. Alternatively, this condition can be expressed 
in terms of the ratio 



which must be the same on both sides of the boundary. It is convenient to refer to 
Z as the "impedance" . The condition that the impedance in both sides of the layers 
interface should be the same gives the possible values for u (eigenvalues). This condition 
is expressed by the equation 



The spectrum in terms of the wave frequency has a continuous part plus a discrete 
part, solution of equation (5.86). The modes w for < u < N 2 are of sinusoidal shape 
in both layers, and for N 2 < uj < Ni the modes w(z) are sinusoidal in the first layer 
and decay exponentially in the second layer. Thus, to deduce how the perturbation 
will change with time from some initial state, it is necessary to represent this state 
as a superposition both of discrete waveguide modes and the continuous spectrum of 
sinusoidal modes. The relative amplitude of the different modes depends on the initial 
state. 



cot 2 (miH) 




(5.86) 



